Parametric resonance of a two-dimensional electron gas under bichromatic irradiation 
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In an ultrahigh mobility 2D electron gas, even a weak nonparabolicity of the electron dispersion, 
by violating Kohn's theorem, can have a drastic effect on dc magnetotransport under ac drive. In 
this paper, we study theoretically the manifestation of this effect in the dc response to the combined 
action of two driving ac-fields (bichromatic irradiation). Compared to the case of monochromatic 
irradiation, which is currently intensively studied both experimentally and theoretically, the presence 
of a second microwave source provides additional insight into the properties of an ac-driven 2D 
electron gas in weak magnetic field. In particular, we find that nonparabolicity, being the simplest 
cause for a violation of Kohn's theorem, gives rise to new qualitative effects specific to bichromatic 
irradiation. Namely, when the frequencies U\ and u>2 are well away from the cyclotron frequency, u> c , 
our simple classical considerations demonstrate that the system becomes unstable with respect to 
fluctuations with frequency i (uii + 1^2). The most favorable condition for this parametric instability 
is I (uii + L02) — lo c . The saturation level of this instability is also determined by the nonparabolicity. 
We also demonstrate that, as an additional effect of nonparabolicity, this parametric instability can 
manifest itself in the dc properties of the system. This happens when ui, L02 and lu c are related as 
3:1:2, respectively. Even for weak detuning between u>i and u>2, the effect of the bichromatic 
irradiation on the dc response in the presence of nonparabolicity can differ dramatically from the 
monochromatic case. In particular, we demonstrate that, beyond a critical intensity of the two fields, 
the equations of motion acquire multistable solutions. As a result, the diagonal rfc-conductivity can 
assume several stable negative values at the same magnetic field. 

PACS numbers: 05.60.-k,73.43.Cd,73.50.Pz 



I. INTRODUCTION 

The cyclotron resonance in a 2D electron gas was first studied almost 30 years agoiS, These studies revealed an 
oscillatory magnetoabsorption of microwave radiation with its principal peak at u) = oj c , where lo c and lo are the 
cyclotron and microwave frequencies, respectively, and several subharmonics at to = nto c due to the disorder-induced 
violation of Kohn's theorem. 

Recently, interest in the properties of microwave-driven electrons in a magnetic field has been revived^ especially 
after experiments^ 5 - carried out on samples with extremely high mobilities, indicated that near the cyclotron resonance 
and its harmonics, irradiation results in drastic changes of the diagonal dc resistivity. In contrast, the Hall resistivity 
remains practically unchanged by illumination, and retains its classical value. The experimental observations^ 5 - were 
confirmed in a number of subsequent gtudies&i&SiifiiAiiiSii&ii This unusual behavior of the weak-field magnetocon- 
ductivity is currently accounted for by an instability resulting from a sign reversal of the diagonal photoconductivity 
under irradiation^*i£ii! The developed instability results in dynamical symmetry breaking, 15 i.e., in an inhomoge- 
neous state of the system characterized by domains of current flowing in opposite directions. Ongoing theoretical 
Studiesi2*i2i22i2ii22i22i2i concentrate on the microscopic description of the sign reversal of photoconductivity. Closely 
related physics was already discussed theoretically quite long agoi^ 5 ^ 

Obviously, a complete understanding of the fascinating properties of ultraclean 2D electron systems under irradiation 
requires additional experimental studies. At the same time, the number of feasible measurements that were not carried 
out so far is limited. A promising avenue seems to be to study the response to bichromatic irradiation. Motivated by 
this, in the present paper we calculate this response within the simple modelSi of a clean classical 2D gas, in which 
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FIG. 1: Dimensionless diagonal conductivity (in units of the Drude value) plotted from Eq. 11411 versus the dimensionless 
magnetic field, defined by Eq. I|15(l . for three cases: (a) monochromatic irradiation with frequency uj\ and dimensionless intensity 
A — 14.4; (b) monochromatic irradiation with the same intensity as in (a) and frequency L02 = 5lui/3; (c) bichromatic case: 
response to simultaneous irradiation with two microwave sources having the same intensities and frequencies as in (a) and (b). 
The emerging region of negative diagonal conductivity is shaded. All three plots (a)-(c) are calculated for | (uj\ + U2) t = 20. 
Full and dashed lines correspond to stable and unstable branches, respectively. 



Kohn's theorem is violated due to nonparabolicity 

i-i- . (i) 



£{p) = 2^ 



2mE 



of the electron dispersion. Here m is the effective mass, and Eq is an energy of the order of the bandgap. 

Denote by £\ and £2 the amplitudes of two linearly polarized ac fields with frequencies uj\ and u)2, respectively. In 
the presence of a dc field Ed Cl the equation of motion for the electron momentum V = p x + ip y takes the form 

& + Z - luJcV + ^-V\V\ 2 = eE dc e* s + ^ (>'* + e~^A + ^ + e~^A , (2) 
at t mEo 2 V / 2 V / 

where lo c is the cyclotron frequency, r is the relaxation time, and 9 is the orientation of the weak dc field with respect 
to the fields £\, £2, which we assume to be parallel to each other. For a monochromatic ac drive, £2 = 0, it was 
demonstrated^ that within a certain interval of magnetic fields near the cyclotron resonance, Eq. J2J yields negative 
diagonal conductivity, <Jd < 0, without significant change of the Hall conductivity. This sign reversal occurs when 
the mobility is high, uj c t 3> 1, and £\ is sufficiently strong. In the vicinity of the cyclotron resonance, <Jd turns 
negative even when the irradiation-induced change of the electron mass is relatively weak. Thus, the simple model 
Eq. @ exhibits negative photoconductivity without invoking Landau quantization. It also predicts bistable hysteretic 
behavior of oa as a function of the detuning from the cyclotron resonance for large enough uj c t. 

In the present paper, we extend the consideration of Ref.|23to the bichromatic case. The most convincing illustration 
that the response to irradiation with two ac-fields cannot simply be reduced to the superposition of the responses 
to each individual field, is presented in Fig. ^ It is seen in Figs. ^a,b) that the individual fields of equal intensity 
and frequency ratio 5 : 3 are unable to reverse the sign of the diagonal conductivity at any magnetic field. At the 
same time, upon simultaneous irradiation by the both fields, a domain of magnetic fields emerges, within which the 
diagonal conductivity is negative (see Fig.^c)). 

In addition, our study reveals the following new features that are specific to the bichromatic case: 

(i) The presence of the second ac field on the r.h.s. of Eq. J2| gives rise to a second domain of magnetic field, within 
which ad is negative. Upon increasing the intensities of the two ac fields, the two domains of negative photoconductivity 
merge into a single domain which broadens much faster with the ac intensity than in the monochromatic case. 

(ii) For monochromatic irradiation, ad could assume either one or two stable values. By contrast, under bichromatic 
irradiation, we find a multistable regime within certain domains of magnetic field. 

(m) In the vicinity of the conditions (lui + L02) = 2lu c and \coi ~ LU2\ — 2uj c , a nonparabolicity-induced parametric 
instability develops in the system. As a result of this instability, the components (u>i + u>2)/2 and \lui — ll>2\/2 emerge 
in addition to the conventional frequencies uj\ and L02 of the momentum oscillations. These components, in turn, 
upon mixing with the components W\, UJ2, give rise to components of V oscillating with frequencies 3a>i — u>2- Thus, 
for bichromatic irradiation, two high-frequency ac driving fields can create a low frequency current circulating in the 
system. In particular, for u>2 = 3wi « 3lu c /2 the system exhibits a dc response to the ac drive. 
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The paper is organized as follows. The case of weak detuning from the cyclotron frequency is treated analytically in 
Sec. [n] In the same section, we present numerical results in this limit, which exhibit nontrivial multistable behavior. 
In Sec. IIIII we consider the case of strong detuning, where we find a nonparabolicity-induced parametric resonance. 
Concluding remarks are presented in Sec. IIVI 



II. WEAK DETUNING 

For monochromatic irradiation, the cyclotron resonance develops when the microwave frequency is close to the 
cyclotron frequency lu c . In this section, we consider bichromatic irradiation when both frequencies uj\ and u 2 are close 
to lj c , \lo\ — lu c \ <C ui c and \u> 2 — ll> c \ -C oj c , so that the cyclotron resonances due to u>i and w 2 can interfere with one 
another. 



A. Calculation of diagonal conductivity 

In analogy to Ref. we search for solutions of Eq. J2J) in the form 

V{t) = Vo + Vi exp(m*ii) + Pf exp(— iw\t) + V 2 exp(iw 2 i) + V 2 exp(— iw 2 i), 



(3) 



where Vq is a small dc component proportional to E Ac . The components V 1 and V 2 are nonresonant and can be 
found from the simplified equations 



~2~' 



(4) 



i(uj 2 + u c )V 2 



e£ 2 
~2~' 



(5) 



where we neglect both relaxation and nonlinearity. However, relaxation and nonlinearity must be taken into account 
when calculating the resonant components and V 2 . Substituting Eq. © into Eq. @, and taking into account 
that \Vi\, \V 2 \ <C \Vi I, \V 2 I , we arrive at a system of coupled equations for the resonant momentum components, 



mEo 



+ |2 



2|P+| 2 ) 



~2~' 



(0) 



t mho 



v? = 



e£i 



(7) 



Despite the inequalities \ <^.\V 1 \ and \V 2 | <C \V 2 |, it is crucial to keep the nonresonant components V 1 and V 2 
when considering the dc component Vq. This yields 



1 2iuj r 

+- + ^{\vt\ 2 + \v. 



t mEo 



2 , ,^ +|2) 



Vo 



2iuo 



mEo 



~ [VtV^ + VtV 2 } V*o = eE A J 



(8) 



Due to the nonlinearity, the microwave intensities induce an effective shift in the resonance frequency lo c . Thus, it is 
convenient to introduce effective detunings f^i and 2 by 



Q 2 = uj 2 - uj c + 



mEo 
mEo 



(|P+| 2 + 2|P+| 



(2\V+ 



and to present formal solutions of Eqs. © and J7J in the form 

e£\T 



2(1 + ifl 1 r) 



Vo 



e£ 2 r 
2 (1 + iQ 2 r) ' 



(9) 
(10) 

(11) 
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Note that the detunings f2i and themselves depend on V± and V^, so that Eqs. 10, (|10|l and i|llfl should be 
considered as a system of nonlinear equations for the resonant momentum components and . Assuming that 
the detunings Sli and fi 2 are known, the solution of Eq. © yields for the dc component 



Po 



eE dc 



(l + tw c T) e ze + 



1 



(efir) (e£ 2 rY 



AmE 1 + iQiT 1 + iQ 2 r 



(12) 



The diagonal conductivity er^ is proportional to Re (Voe~' ie ). Thus, the second term in Eq. (|T^|) gives rise to a 9- 
dependence of the nonparabolicity- induced contribution to the diagonal conductivity which is given by sin (29 — <p). 
Here, cf> satisfies the equation 



tan 6 = 



£7it(1 



mr 2 ) 



E\ + n 2 r (1 + tt 2 r 2 ) El 



i + n 2 2 r 2 ) £\ + (i 



(13) 



We deduce that the minimal value of Ud is given by 



ne 



AmEn 



(£fn 2 T + El^rf + (El + El) 
(1 + ft 2 r 2 ) (1 + ft 2 r 2 ) 



1/2. 



(14) 



In the following, we will be particularly interested in <r 
of the zero-resistance state. 



since the condition a d nm < is sufficient for the formation 



B. Numerical results: Multistability 



As demonstrated in Ref. 1271 the diagonal conductivity in the monochromatic case shows a region of bistability. 
In the bichromatic case under study, there can even be multistable behavior as will now be shown. We measure the 
frequency difference of the ac- fields by A = (u)\ — lo 2 )t and the magnetic field by 



b = 



LJl + U) 2 



(15) 



which depends linearly on the magnetic field B. Upon substituting the formal solutions and of Eq. (|llfl into 
Eqs. (|5ITX)|) . these can be written as a pair of coupled equations for the effective detunings f2ir and fl 2 r 



1 2 



,4 



2t7 2 



i + ^tY i + (n 2T y 



f2 2 r = b + A 



i + (n lT y i + (n 2T y 



where 77 and A, given by 



77 = £ 2 /£i 



A = LU C T 



(e£irf 



AmEi 







(16) 



(17) 



measure the ratio of the field amplitudes and the ratio of the absolute field intensities to the nonparabolicity of 
the electron spectrum, respectively. As in the monochromatic case, the strength of the first order correction to the 
Drudc conductivity is proportional to the microwave intensity and thus A. At a fixed magnetic field b and at fixed 
ac frequencies and amplitudes, this coupled system of two third-order equations can yield up to nine simultaneous 
solutions (f^iT, f^i") for the effective detunings. Since a d nm is directly related to these effective detunings via 



<JD { 1 



A 

LO c T 



1 



2rf (1 + QiTfl 2 T) 



1 + (r>!T) 2 1 + (n 2 r) 2 (1 + ( r> lT ) 2 ) (i + ( r> 2r )- 



1/2' 



> - 



(18) 



where ct_d = ne 2 /(oj 2 t) is the Drude conductivity, there are thus up to nine individual branches of a™ m at a given 
b. This multistable behavior occurs in the vicinity of the resonance and will be studied below for the specific case of 
77 = 1. 

We first focus on the dependence of cr d nm (b) on A. For large A, i.e. markedly different ac-frequencies, we expect 
two separate regions in b where <j™ m deviates significantly from the Drude result. These are the regions where the 
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b b b 

(a) (b) (c) 

FIG. 2: Evolution of the dimensionless (in units of the Drude conductivity cm) minimal conductivity <r™ m as function of 
magnetic field b, defined in Eq. 1151k for three different values of A; (a) A = 20, (b) A = 5, (c) A = 1. The curves are 
calculated for the values of parameters r] = 1, A = 5 (defined by Eq. 1171 '). and u c r = 25. The distance of the two dips that 
can be clearly discerned in (a) is roughly A. When lowering A, the dips move closer together (b) and finally merge (c). In 
addition, multistable regions emerge. As in Fig. unstable branches are plotted as dashed lines. 



cyclotron frequency is in resonance with one of the two ac frequencies, i.e. either u)\ ~ oj c or 0J2 — to c . Inside these 
regions, the behavior with respect to b is very similar to the monochromatic case, except that the irradiation-induced 
effective shift of u c now depends on both external frequencies. In particular, the emergence of bistable regions inside 
these two separate intervals as in the monochromatic case is to be expected. This can be seen in Fig. |2Ja), where 
a mm j g ghown as a function of magnetic field b for rather large A. Two dips in <j™ m can be clearly discerned, the 
inner branches of which are unstable. Upon reducing A, the two dips move closer together up to a point where the 
frequencies uj\ and UJ2 are so close that the analogy to the monochromatic case breaks down and the two dips start 
to interact to finally from a single multistable dip in the limit A — > 0. This behavior is exemplified in Figs.|2fb,c). It 
can be seen that multiple solutions of <7™ m (&) develop upon reducing A. 

Next, we consider the case of negative diagonal conductivity and study the evolution of <j™ m with magnetic field. 
As expected, there is a threshold intensity below which no negative diagonal conductivity is observed. Upon increasing 
the field amplitudes and thus A, the negative first order correction to the Drude conductivity grows linearly with 
A/{lo c t) as can be seen from Eq. I|18(l . When this correction exceeds one, negative cr™ m is to be expected in some 
regions of magnetic field b. Fig. [3] shows the 6-dependence of cr™ m for three specific values of A. It can be seen that 
at low A, no regions in b with negative diagonal conductivity can be observed. At higher A, two regions in b show 
negative cr™ n -branches as is also indicated by the shaded regions in Fig.|3Jb). For even higher A, a single large region 
in b shows negative diagonal conductivity. 

To clarify the evolution of these regions with increasing field amplitudes, we plotted the extension of the regions in 
b as a function of A. The result is shown in Fig. It is remarkable that above the threshold value of A, first a single 
region appears that shows negative a™ m . Then, in the immediate vicinity of the threshold a second, well separated 




FIG. 3: The dimensionless minimal diagonal conductivity, a™ m (b), denned by Eq. 1181 . is plotted for three different values of 
the dimensionless intensity A of the two ac-fields; (a) A — 20, (b) A = 30, (c) A = 32. The domains of magnetic field b with 
negative ctJ"™ are shaded. The dotted line indicates the boundary between positive and negative cr™ n . Unstable branches are 
dashed as in Fig. The curves are calculated for the values of the parameters r) = 1, A = (uji — lu2)t = 25 and ui c r = 25. 
In (b), the point P is a point where a continuous stable and a continuous unstable branch intersect without "noticing" each 
other, while the point P' is a cusp which separates a stable from an unstable branch. 
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FIG. 4: Evolution of the regions of negative <rj" n with irradiation intensity, A. Shown are the Dichromatic case (left panel) 
and the monochromatic case (right panel). In both cases, only the dominant (stable) branches are shown to avoid confusion. 
The parameters are the same as in Fig. [3] 



region develops. Upon further increasing A, the width of these regions grows and, eventually, the two regions merge 
to form a single broad region of negative diagonal conductivity. For comparison, we also show the monochromatic 
case in the right hand panel of Fig. 0] 



C. Stability of different branches 



The stability of the various branches of <j™ m (b) as shown by solid and dashed lines in Figs. 11121 and|3 can be 
obtained from a standard stability analysis. As usual, we find that stable and unstable branches "meet" at cusps, 
as clearly seen at the minima of a a in Figs. ^ and [21 The transitions from unstable to stable branches at larger b in 
these figures are also accompanied by cusps, although this can not necessarily be discerned within the resolution of 
the figures. The number of branches increases with the irradiation intensity, cf. Figs. |21 and |21 The rule that stable 
and unstable branches meet in cusps remains valid, although this statement becomes less trivial. For example, in Fig. 
E^b) stable and unstable branches intersect at point P without "noticing each other" . Accordingly, there is no cusp 
at this point. At the same time, there is a cusp at P' in Fig. E^b) where the same branches switch between stable 
and unstable. Figs. |2Ic) andOUc) illustrate how new branches and multistability emerge with increasing irradiation 
intensity. The emergence of new stable and unstable branches occurs in pairs which meet at additional cusps. In 
both figures |2Ic) and^c), there are regions in magnetic field with three coexisting stable solutions (tristability). 
Further increase of A would lead to up to eight cusps in Fig. OJc), each of which is a meeting point of stable and 
unstable branches. Thus, the tristability situation illustrated in Fig. |2{c) will evolve into a magnetic field domain 
with "four-stability" . 



III. STRONG DETUNING 



In this section, we consider the case when both frequencies lo\ and cj 2 are tuned away from uj c . This implies that 
the system (|6I7|) decouples and acquires the obvious solutions 

-p+ = i j}+ — i (19) 

1 2 [1 + i (wi - cj c ) r] ' 2 2 [1 + i (cu 2 -lo c )t}' v ' 

The condition e£±, e£i <C to c (mEy) 1 ^ 2 for decoupling follows from Eqs. (|6I7[) . assuming that \u>x — u) c \, \uj2 — w c | 
~ tu c ^> 1/r. Interestingly, even under this condition, the solutions are unstable for certain relations between the 
frequencies u)\, 102- The mechanism for this instability relies on mixing of the two external drive frequencies by the 
nonparabolicity which results in a modulation of the effective cyclotron frequency. This modulation, in turn, can lead 
to 'parametric resonance. 

To perform the stability analysis of the solutions Eq. (|19f) . we introduce a small deviation V — > P + SV and linearize 
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Eq. J5J with respect to 8V. The linearized equation (J2J has the form 

d 



dt 



(SV) + - - iu c + ^|7>| 2 <5P + -^7> 2 = 0. 



(20) 



This equation couples SV to (SV) via the nonparabolicity of the electron spectrum. The corresponding equation for 
SV* reads 



dt 



(SV*) 




^\ v] A S v*-^r(V*) 2 SV = 0. 
mEn / mhn 



(21) 



The coupling coefficient, V 2 , as seen from Eq. @, contains the harmonics ±2oji, ±2^2, ±(^i +W2), and ±(wi — u^)- 
This suggests that SV(t) also contains a number of harmonics, namely, ±wi, ±w 2 , ±(^i + ^2)/2, and ±(wi — w 2 )/2. 
An instability might develop when one of these frequencies is close to w c . Thus, in the monochromatic case, the 
instability develops only in the vicinity of the cyclotron resonance uj\ ss oj c . The branches, shown with dashed 
lines, in Figs. ^a,b), are unstable due to this instability. By contrast, the bichromatic case offers two additional 
options for an instability to develop, even if the frequencies u)±, 102 are nonresonant, namely uj c ~ (^1 + ^2) /2 and 
uj c w I (oji — LJ2) |/2. The considerations of both cases are analogous to each other. Therefore, we focus on the first 
case below. 



A. Parametric instability at (cui +0^2) ~ 2w c 



Upon substituting the ansatz 



SV(t) = Cexp 



i (loi + oj 2 ) 



SV(t)* = C* exp 



i (uji + uj 2 ) 



(22) 



into Eqs. I|2UI21|I and keeping only resonant terms, we obtain the following system of algebraic equations for C and 

C* 



r + i 

r 

T 



i (u>i -\- LO2 — 2lu c ) 2ilo c 



i (u>i + UJ2 — 2ui c ) 



mEo 
2iuj c 



(\n\ 2 
wt\ 2 ^ 



c = 



wt\ 2 ) 



c* 



2iu c 
2iu. 



vtv^c* 



mEn 



{vtvtyc. 



2 mE 

Thus, the most favorable condition for instability is determined by the following relation between oj\ and L02 

2 



(23) 
(24) 



2ui r 



2uj r 



mEo 

,,2 



(]vt\ 2 + wt 



2mE n 



f 2 



(UJI - CJ C ) (UJ2 ~ UJ C ) 



2uOr 



2e 2 {£l + £l) 
mEo (wi — UJ2Y 



where Eq. 119|) has been used. In this case, the increment T is maximal and given by 



l + (*H!L\ \V + V + \^- 1 - 



e 2 ui c £\£2 



mEo (Ull - LU C ) (UJ2 - LUc) 



e 2 (loi + lu 2 ) \E\£-2 
mEo (uj\ - uj 2 ) 2 



(25) 



(26) 



The parametric instability develops if T rnax is positive. It is important to note that the condition T max > is 
consistent with the condition of strong detuning when the simplified expressions Eq. are valid. Indeed, assuming 

|cJi — uj c \ ~ \uj2 — w c | ~ oj c , the two conditions can be presented as ui c (mE ) 1/2 > e£i,e£ 2 > T" 1 (mE ) 1/2 . Therefore, 
for lu c t 3> 1, there exists an interval of the amplitudes of the ac fields within which both conditions are met. Note 
also, that parametric resonance does not develop exactly at uj c — (uj\ + lu 2 )/2, i.e. at b = (in dimensionless units, see 
Eq. In fact, from Eqs. I|25|l and (|26|l it can be concluded that T max > corresponds to b > 1. In experimental 

situations, when loi and 102 are fixed, Eq. (|25|l can also be viewed as an expression for the magnetic field uj c = uj° pt , 
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at which the parametric instability is most pronounced. The interval of lo c around u° pt , within which the increment 
is positive can be found from the dependence T (cj c ) 




opt 



\ 



e 2 (oji + lj 2 ) \£i£ 2 \ 



mE (lui - lu 2 ) 

Upon setting T(uj c ) — in the l.h.s. of Eq. (|27[) . we find the width of the interval to be 



(27) 



\u c -^\ < 



e 2 Lo c \£i£ 2 \ 



2mEn (lui — uj c ) 



T '1 


r-i 







e 2 (wi + lo 2 ) \£i£ 2 \ 
mE (ll>i - ui 2 ) 2 



1/2 



(28) 



It is instructive to reformulate the condition for the parametric instability in a different way. Assume for simplicity 



that £\ = £ 2 . Then the combination e \£i£ 2 \/2mEn(uJi 



is equal to 5m/m, where 5m /m is the relative 



correction to the electron effective mass due to irradiation. 27 From Eq. (|26|l it follows that the condition T max > 
can be presented as (uj c t) (Sm/m) > 1. With lu c t 1 this condition can be satisfied even for 5m <C m. Note, that 
in the case of weak detuning, the same condition is required for the dc conductivity to assume negative values near 
the cyclotron resonance. 

Summarizing, we arrive at the following scenario. In the case of strong detuning, there is no mutual influence of 
the responses to the ac fields £\ and £ 2 as long as they are weak. However, as the product \£\£ 2 \ increases and 
reaches a critical value \£\£ 2 \ c ^ the threshold, T max = 0, where T max is given by Eq. I|2t)|) . is exceeded at the magnetic 
field uj c = ui° pt determined by Eq. (|2"5|) . Above the threshold, fluctuations with frequencies close to (oji + uo 2 ) /2 are 
amplified. This effect of parametric instability is solely due to nonparabolicity. Then the natural question arises: At 
what level does the component of momentum with frequency {oj\ + uj 2 ) /2 saturate above threshold? This question is 
addressed in the next subsection. 



B. Parametric instability at \uo\ — U2\ ~ 2uj c 

We now briefly discuss parametric instability at weak magnetic field, lo c « | (ui — oj 2 ) /2|. Assume for concreteness, 
that uj\ > ijj 2 . In this case, the optimal magnetic field, Cj° pt , is lower and reads 



UJl - UJ 2 



2Cj° pt 



u 



£; 



£l 



(wi - Cj° pt ) 



(ljj 2 - uj° Pt ) 




2e 2 
mEn 



£? 



£l 



(lji+uj 2 ) 2 {2>uj 2 -u>i) 2 



At 



j° pt , the threshold condition for parametric instability, analogous to Eq. (J2HJ), has the form 



1 e 2 (uji - u> 2 ) \£i£ 2 \ 

t mEo (u>i + u)2) |3a;2 — wi| 



> 0. 



(29) 
(30) 



There is no real divergence in Eqs. (|29|l . (|30|l in the limit uj\ — » Zlo 2 , since these are derived under the assumption 
that the difference \2>u 2 — u\ \ is > 1/r. 



C. Saturation of parametric resonance 



As the threshold for parametric resonance is exceeded, the harmonics with frequency (wi + lo 2 ) /2 can no longer be 
considered as a perturbation, but rather have to be included into the equation of motion. In other words, we must 
search for a solution of Eq. @ in the form 



V = Vt exp (iuit) + Vt exp (iu 2 t) + V 3 (t) exp 



(31) 



where Vs(t) is a slowly varying function of time. Upon substituting this form into Eq. J5J), we obtain the following 
coupled equations for Vz{t) and V^(t) 



dV 3 
dt 



1 i (lo\ + uj 2 — 2lu c ) 
t + 2 



mEn 



+ |2 



217^. 



+ |2 



\r,\ 2 ) 



2lLU c 

mEc 



'V + 'P + V* 
1 1 '2 1 3 i 



(32) 
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dt 



i (u>i + u>2 — 2lo c ) 



mE Q 



{2\r+\ 2 + 2\r+\ 2 + \v 3 \ 2 ) 



6 mE y 1 2 



V*. 



(33) 



Saturated parametric instability is described by setting dV^/dt = and dV^/dt = in Eqs. and (JJSJl, respectively. 
The result for V$ has the simplest form for the optimal magnetic field uj c 



Opt 



1^3 1 («. 



opt \ _ 



\V+P+\ 2 - 



4m 2 Eq 



1/4 



(Wl +0J 2 ) T 2 



From Eq. (|34[1 we conclude that, in the vicinity of the threshold, \Vs\ increases as ( | £T X £" 2 1 — \£i£ 2 



1/4 



(34) 



oc r 1/4 



Well above the threshold it approaches the value y\V^P 2 \. From here we conclude that, even upon saturation, the 
magnitude of the nonparabolicity-induced harmonics with frequency (lui + lo 2 ) /2 does not have a "back" effect on 
the magnitudes Eq. I|19|l of the responses to the oc fields. 

For magnetic fields in the vicinity of uj° pt the saturation value, |"p3|(a; c ), is given by 



1^3 IK) 



4m 2 E 2 



(lui + uj 2 ) t 2 



1/2 



2mEo\Lu c 



opt 1 



1/2 



(lOi + LJ 2 ) 



|P3| 2 K° Pt ) 



2mEo\ui c — uj\ 



opt \ 



(u)i + UJ 2 ) 



1/2 



(35) 



/ \!/2 

In contrast to \P 3 \(uj° pt ) , the threshold behavior of I^Ku^) is slower, namely I^K^c) oc (|fi^2| — |f 1^2)0 ) • In 

principle, one has to verify that the solutions Eqs. lf3*5|l . that describe the saturated parametric resonance, are 

stable. This can be done with the use of the system Eqs. by perturbing it around the saturated solution. 

The outcome of this consideration is that the corresponding perturbations do indeed decay. 



D. Implications for dc transport 

As we demonstrated in the previous subsection, the parametric instability that develops in the case of two ac fields 
above a certain threshold, results in the component of the momentum, P3 exp [i (uji + L02) £/2] , which well above the 
threshold has the magnitude 7^3 ~ \/\P\\\P2\- The important consequence of the developed parametric resonance 
is, that the component Vz gives rise to new harmonics in the term oc \V\ 2 V in the equation of motion Eq. (J2J. Of 
particular interest are the T-Vinduced terms 



\P\Vt + V1VZV3] exp 



. 3wi 



■ w 2 



(36) 



It is easy to see that, under the condition lu 2 ~ 3wi, these terms act as an effective dc field, and thus generate 
low-frequency circular current even without dc drive. If the relation between the frequencies is precisely 1 : 3, then 
the magnetic field, at which the developed parametric instability would give rise to a quasistationary circular current 
distribution, can be determined from Eq. (f 251) 



uJr 



2wi 



2e 2 (£ 2 



£1) 



mE (wi — uj 2 y 



(37) 



If the ratio uo 2 /oji is close, but not exactly 1 : 3, there is still a certain allowance, determined by Eq. I|28|l for 
the formation of the quasistationary current. The above effect of spontaneous formation of dc-like currents under 
irradiation is distinctively different from the formation of current domains when aa turns negative under irradiation. 
Firstly, the effect is specific to bichromatic irradiation. Secondly, it requires rather strict commensurability between 
the two frequencies, and finally, it develops within a very narrow interval around a certain magnetic field. 



IV. CONCLUSIONS 



In the present paper we have considered the problem of single electron motion in a magnetic field under irradiation 
by two monochromatic fields. When the frequencies lo± and lj 2 differ only slightly, \u>\ — lo 2 \ ~ 1/r, the effect of a weak 



10 



nonparabolicity of the electron spectrum on the diagonal conductivity is qualitatively the same for monochromatic^ 
and bichromatic irradiation. The prime qualitative effect which distinguishes the bichromatic case is the emergence 
of a parametric resonance at magnetic fields uj c = (uji + u>2)/2 and lo c — \u)\ — u>2\/2 when the detuning is strong 
(of the order of the cyclotron frequency). It is instructive to compare this effect with the parametric resonance of 
electrons in a magnetic field due to a weak time modulation of the field amplitude i^&LQ The latter effect, considered 
more than 20 years ago, has a transparent explanation. The modulation of the magnitude of a dc field with frequency 
2uj c translates into a corresponding modulation of the cyclotron frequency, so that the equation of motion of the 
electron reduces to that for a harmonic oscillator with a weakly time-modulated eigenfrequency. The solution of 
this equation is unstable, if the modulation frequency is close to 2uj c . As a natural stabilizing mechanism of the 
parametric resonance, the authors of Ref. considered the nonparabolicity Eq. of the electron dispersion. For a 
characteristic magnetic field of B = 0.1 T, the cyclotron frequency is uj c — 3.6 • 10 11 Hz, i.e. in the microwave range 
so that conventional modulation of B with a frequency 2lo c is technically impossible. To bypass this obstacle, it was 
proposed in Ref. to use microwave illumination with frequency 2u> c to create a parametric resonance. The idea 
was that the magnetic field component of the pumping electromagnetic wave would provide the necessary oscillatory 
correction to the external dc magnetic field. In the present paper, we ha ve dem onstrated that two nonresonant ac 
sources can enforce a parametric resonance of the type considered in Rcfs. 29 30 31 without any time modulation of 
the dc magnetic field. Remarkably, this bichromatic-radiation-induced cycl otron reso nance emerges due to the same 
nonparabolicity Eq. that played the role o f a stabilizing factor in Refs. I29l30l3ll Roughly, the time modulation 
of uj c in the dc field required in Refs. t29l30l3ll for parametric resonance emerges from the "beatings" of the responses 
to the two ac signals. The nonparabolicity transforms these beatings into a modulation of the cyclotron frequency. 
Although the increment, T, for parametric resonance, induced by bichromatic microwave irradiation, is proportional 
to the product £\£2 of the amplitudes of the two sources, while in Ref. OH] it was proportional to the first power of 
the magnetic component of the ac field, the "bichromatic" increment is much bigger. As demonstrated above, the 
bichromatic increment is T ~ lj c (mc 2 / Eq) [£±£2/ B 2 ), which should be compared to the increment r ~ uj c (£/B) of 
Ref. l3ll The ratio contains a small factor (£/B) which is offset by the huge factor (mc 2 /Eq). 
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